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The nonequilibrium dynamic phase transition, in the two 
dimensional site diluted kinetic Ising model in presence of an 
oscillating magnetic field, has been studied by Monte Carlo 
simulation. The projections of dynamical phase boundary 
surface are drawn in the planes formed by the dilution and 
field amplitude and the plane formed by temperature and field 
amplitude. The tricritical behaviour is found to be absent in 
this case. 

PACS number(s): 05.50.+q 



I. INTRODUCTION 

Though the Ising model was proposed nearly three 
quarters of a century ago its dynamical aspects are still 
under active investigation. Nowadays, the study of the 
dynamics of Ising models in presence of time varying 
magneic field, became an active and interesting area of 
modern research. The dynamical response of the Ising 
system in presence of an oscillating magnetic field has 
been studied extensively by computer simulation |^-^] in 
the last few years. The dynamical hysteretic response 
[ 1 - 3 1 and the nonequilibruim dynamical phase transition 
[ 4 — S | are two important aspects of the dynamic response 
of the kinetic Ising model in presence of an oscillating 
magnetic field. 

Tome and Oliviera [Q first studied the dynamic tran- 
sition in the kinetic Ising model in presence of a sinu- 
soidally oscillating magnetic field. They solved the mean 
field (MF) dynamic equation of motion (for the aver- 
age magnetisation) of the kinetic Ising model in presence 
of a sinusoidally oscillating magnetic field. By defin- 
ing the order parameter as the time averaged magneti- 
sation over a full cycle of the oscillating magnetic field 
they showed that the order parameter vanishes depend- 
ing upon the value of the temperature and the amplitude 
of the oscillating field. Precisely, in the field amplitude 
and temperature plane they have drawn a phase bound- 
ary separating dynamic ordered (nonzero value of order 
parameter) and disordered (order parameter vanishes) 
phase. They B have also observed and located a tri- 
critical point (TCP), (separating the nature (discontinu- 
ous/continuous) of the transition) on the phase boundary 
line. 

Since, this transition exists even in the static (zero fre- 
quency) limit such a transition, observed B from the so- 



lution of mean field dynamical equation, is not dynamic 
in true sense. This is because, for the field amplitude less 
than the coercive field (at temperature less than the tran- 
sition temperature without any field), the response mag- 
netisation varies periodically but asymmetrically even in 
the zero frequency limit; the system remains locked to 
one well of the free energy and cannot go to the other 
one, in the absence of noise or fluctuation. 

Lo and Pelcovits || first attempted to study the dy- 
namic nature of this phase transition (incorporating the 
effect of fluctuation) in the kinetic Ising model by Monte 
Carlo (MC) simulation. In this case, the transition dis- 
appears in the zero frequency limit; due to the presence 
of fluctuations, the magnetisation flips to the direction 
of the magnetic field and the dynamic order parame- 
ter (time averaged magnetisation) vanishes. However, 
they H have not reported any precise phase boundary. 
Acharyya and Chakrabarti |6[] studied the nonequilib- 
rium dynamic phase transition in the kinetic Ising model 
in presence of oscillating magnetic field by extensive MC 
simulation. It was also noticed by them |1 that this dy- 
namic phase transition is associated with the breaking of 
the symmetry of the dynamic hysteresis (m — h) loop. 
In the dynamically disordered (value of order parame- 
ter vanishes) phase the corresponding hysteresis loop is 
symmetric, and loses its symmetry in the ordered phase 
(giving nonzero value of dynamic order parameter). They 
have Halso studied the temperature variation of the ac 
susceptibility components near the dynamic transition 
point. They observed that the imaginary (real) part 
of the ac susceptibility gives a peak (dip) near the dy- 
namic transition point (where the dynamic order param- 
eter vanishes). The conclusions were: (i) this is a distinct 
signal of phase transition and (ii) this is an indication of 
the thermodynamic nature of the phase transition. 

It may be mentioned here that the statistical distri- 
bution of dynamic order parameter has been studied by 
Sides et al Q. The nature of the distribution changes 
near the dynamic transition point. They have also ob- 
served Q that the fluctuation of the hysteresis loop area 
becomes considerably large near the dynamic transition 
point. 

Very recently, the relaxation behaviour, of the dy- 
namic order parameter, near the transition point, has 
been studied |§] by MC simulation and solving meanfield 
dynamic equation. It has been observed that the relax- 
ation is Debye type and the relaxation time diverges near 
the transition point. The 'specific heat' and the 'suscep- 
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tibility' also diverge |9|] near the transition point in a sim- 
ilar manner with that of fluctuations of order parameter 
and fluctuation of energy respectively. 

In this paper, the dynamic phase transition has been 
studied in the site diluted (by nonmagnetic impurities) 
kinetic Ising model by MC simulation. The phase bound- 
aries are plotted in the planes formed by the field ampli- 
tude and the temperature and in the plane formed by 
the impurity concentration and the field amplitude. The 
paper is organised as follows: in section II the model and 
the simulation scheme are discussed, the MC results are 
given in section III and the paper ends with a summary 
of the work in section IV. 



II. THE MODEL AND THE SIMULATION 
SCHEME 

The Hamiltonian, of a site diluted Ising model (with 
ferromagnetic nearest neighbour interaction) in presence 
of a time varying magnetic field, can be written as 

H = - Jiisfs'cjCj - h(t) £ sfci. (2.1) 

<ij> i 

Here, sf(= ±1) is Ising spin variable, Jy is the inter- 
action strength, a (— or I) represents the site (i) 
which is either occupied (cj = 1) or vacant (ci = 0). 
h(t) = hocos(uit) is the oscillating magnetic field, where 
ho and to are the amplitude and the frequency respec- 
tively of the oscillating field. The system is in contact 
with an isothermal heat bath at temperature T. For 
simplicity all Jij and the value of the Boltzmann con- 
stant are taken equal to unity. The boundary condition 
is periodic. 

A square lattice of linear size L(= 100) has been con- 
sidered. The lattice sites are randomly occupied by mag- 
netic sites with a finite probability p. So, the degree of 
dilution or the concentration of (nonmagnetic) impuri- 
ties, is q = 1 — p. At any finite temperature T and for a 
fixed frequency (uj) and amplitude (ho) of the field, the 
dynamics of this system has been studied here by Monte 
Carlo simulation using Metropolis single spin-flip dynam- 
ics. Each lattice site is updated here sequentially and one 
such full scan over the lattice is defined as the time unit 
(Monte Carlo step per spin or MCS) here. The initial 
configuration has been chosen such that the all spins are 
directed upward. The instanteneous magnetisation (per 
site), m(t) = (1/L 2 ) ^ s\ci has been calculated. From 
the instanteneous magnetisation, the dynamic order pa- 
rameter Q = § m(t)dt (time averaged magnetisation 
over a full cycle of the oscillating field) is calculated. This 
dynamic order parameter is a function of temperature 
(T), field amplitude (ho) and the impurity concentration 
(q) ,i.e, Q = Q(T,ho,q). Each value of Q has been cal- 
culated by averaging over 25 number of initial impurity 
realisations. The frequency of the oscillating magnetic 
field used here is equal to 0.0628. 



III. RESULTS 

It has been observed that Q = Q(T,ho,q) is nonzero 
for a finite set of values of ho,T and q, and Q vanishes 
elsewhere. In the space formed by h ,T and q, there is 
a surface which divides the Q — region from Q / 
region. Fig. I shows the schematic diagram of this phase 
surface. 

Previously, a number of numerical studies [4,6,8] are 
performed regarding the dynamic transition in the q = 
plane (i.e., projection of this phase surface on ho — T 
plane. In that case, it was observed that in the ho — T 
plane there is a distinct phase boundary below which 
Q is nonzero and above which Q vanishes. There is a 
tricritical point on the phase boundary which separates 
the nature (discontinuous/continuous) of this transition. 

Fig. 2a shows the dynamic phase boundary in the ho — 
T plane for different values of the impurity concentration. 
It has been observed that as the impurity concentration 
increases the phase boundary shrinks inward. In this 
case, the entire phase boundary has been scanned and 
the transition observed is always continuous. Unlike the 
earlier case Q , no such tricritical point is observed here. 
Fig. 2b shows the temperature variations of the dynamic 
order parameter (i.e, Q versus T) for two different values 
of ho in the case of very weak disorder (inpurity). 

A similar kind of dynamic phase boundary has been 
obtained in the ho — q plane (i.e., the projection of the 
phase surface on ho — q plane). Fig 3a. shows the phase 
boundaries in the hg — q plane for different values of tem- 
peratures. Like the earlier case, here also, as the tem- 
perature increases the phase boundary shrinks inward. 
Here also the transition is continuous along the entire 
phase boundary. Two typical transitions (fall of order 
parameter with respect to the impurity concentrations) 
are shown in Fig. 3b. 

IV. SUMMARY 

The nonequilibrium dynamic phase transition, in the 
site diluted kinetic Ising model in presence of oscillating 
magnetic field, is studied by Monte Carlo simulation. 

The value of the dynamic order parameter gets nonzero 
below a boundary surface in T, ho and q space, and above 
the surface it vanishes. The projections of this surface on 
ho—T plane and on ho — q plane are drawn. The nature of 
the transition observed here is continuous along the en- 
tire phase boundary in the ho — T plane with very small 
impurity concentrations. This is unlike the case observed 
earlier || in the pure sample. A similar kind of transition 
is observed for a fixed temperature with varying impu- 
rity concentrations. Here also no tricritical behaviour 
was observed. The studies reported here are mostly ob- 
servational, no attempt has been made to understand 
these phenomena from the knowledge of the theoretical 
background. 
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There is some experimental evidence of dynamic tran- 
sition. Recently, Jiang et al [jlo| observed the dynamic 
transition in the ultrathin Co/Cu(001) sample by mag- 
neto optic Kerr effect. They have noticed the dynamical 
symmetry breaking of the hysteresis loop. However, the 
detailed quantitative study to draw the phase boundary 
is not yet done. 

It is quite interesting to study the dynamic transition 
in randomly varying (in time but uniform in space) mag- 
netic field. The work in this front is going on and 
the results will be published elsewhere. 



parameter (Q) for two different values of field amplitudes 
(ho). 

Fig. 3. (a) Projections of dynamic phase surface on the 
ho — q plane. (O) represents T — 0.25, (+) represents T 
= 0.50 and (□) represents T = 0.75. (b) Variations of 
dynamic order parameter (Q) with respect to impurity 
concentration (q) for two different values of field ampli- 
tudes (ho). 
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Figure captions 



Fig. 1. Schematic diagram of dynamic phase boundary 
surface in the space formed by T, ho and q. Below this 
surface Q ^ and above the surface Q = 0. 

Fig. 2. (a) Projections of dynamic phase surface on the 
h — T plane. (O) represents q = 0.05 and (+) represents 
q = 0.3. (b) Temperature variations of dynamic order 
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